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(— | ' Abstract. Supersymmetry, shape invariance, exact solubility, and the factorization 

method are often studied together in the literature. At the dawn of these topics 
confusion was present in regards to their scope of applicability and the relation among 
them. Considerable work have been put to study and resolve the relation among two 
or more of these topics. These works are scattered over the literature. While looking at 
the literature, one can not overlook the number of places where authors confuse these 
terms, and concluding implications depending on wrong assumptions of the relation 
between two or more of these topics. In this letter we define supersymmetry, and 
\q ' shape invariance, and show the relations which connects them to exact solubility and 

, the factorization method, referring to the literature for the respective detailed work 

C*^ ' and proofs. At last we conclude our letter with a Venn diagram which illustrates those 

relations. 



1. Supersymmetric Quantum Mechanics 

Witten, in his seminal paper [I], defined the algebra of a supersymmetric quantum 
system, which he derived from the algebra of supersymmetry in field theory. In 
supersymmetric quantum systems, there are charge operators Qi which commute with 
the Hamiltonian 

[Qi,H)=0, i = l,2,...,JV (1) 

and they obey the algebra 

{Qi,Qj} = ^H (2) 

where TC is called the Supersymmetric Hamiltonian. Witten stated that the simplest 
quantum mechanical system has N = 2, it was later shown that the case where N = 1, 
if it is supersymmetric, it is equivalent to an = 2 supersymmetric quantum system 
[2]. In the case where iV = 2 we can define 

Q = (-L)(Q 1+i Q 2 ) 7 Qt = ( _L )(gi _,g 2) (3) 
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these operators will obey the algebra 

h = {q,q% g 2 = g t2 = o (4) 

consequently, 

[Q,H] = [Q^H] = (5) 
We can represent our charge operators as 2 x 2 matrices 

Q=(°l) M 
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where A is a linear operator, and A* is its adjoint. Now the supersymmetric 
Hamiltonian has the representation 

where in the last step we designated the hermitian operators A^A and AA' as the 
Hamiltonians H_ and H + , respectively. It is important to notice that H_ and H + are 
positive semi-definite, for if \a ) is an eigenstate in the Hilbert space of any of these 
Hamiltonians, say H_ then 

(a\H_\a) = a (a\a) = a (8) 

where a is the eigenvalue of \a ). At the same time 

(a\H-\a) = (a\A^A\a) = \A\a)\ 2 (9) 

where |.| is the norm of our Hilbert space, whose range is always real. Since the 
right hand side of (IE]) must be equal to the right hand side of (J9j), and by the positivity 
axiom of our metric space (Hilbert space), \\x )| > for all the vectors in our metric 
space [3l pp.14], so it follows that if_ is positive semi-definite. It also follows that H 
has eigenvalue zero iff A annihilates a vector, which must be the ground state. The 
same argument can be applied to H + . For an alternative proof see [H Th. X25, pp 180]. 

It is straight forward to show that if ip~ is an eigenstate for H^ with a positive 
eigenvalue E~ 

H_^p- = A^A^/j- = E~ip- (10) 

then tp + = (E~)~ l l 2 (Ail)~) is a normalized eigenstate for H + with an eigenvalue 
E+ = E~ 

H + iP + = —L=AAiAip- = E-iP + (11) 
v E 
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in the same manner transforms an eigenstate of H + to an eigenstate of ff_ with 
the same energy [5], [6]. This can be written as 

(12a) 




(12b) 

This result comes with no surprise for us, since we constructed our charge 
operators (CD) to obey our algebra (jlj), and we saw that this algebra implies that the 
charges commute with the Hamiltonian ([5]), this commutativity is responsible for this 
degeneracy. 

In literature, the linear operators A and A' are usually considered in the following 

form 

A = 4- + w(x) (13a) 
dx 

A* = — f- + w(x) (13b) 
dx 

where w(x) is an arbitrary function of x and it is called the superpotential. Now, 
we can look at the differential form of our partner Hamiltonians 

H_ = A ] A = —4— + w 2 (x) - w'(x) (14a) 
dx 2 

H + = AA ] = -4-z + w 2 (x) + w'(x) (Ub) 
dx 2 

the prime is dj dx, and we set fi = 2m = 1. The potentials in the above Hamiltonians 
are called the supersymmetric partner potentials or superpartners, and they are given 
by 

V± = w 2 (x) ±w'(x) (15) 

this nonlinear differential equation is known as Ricatti equation. A node-less 
solution to the Hamiltonians can written in terms of the superpotential 



i)^ = Nexp(±J w(t)dtj 



(16) 



where ipQ is a solution of H_ and ipQ is a solution of H + , the subscript is the 
number of nodes in the function, zero in this case. iV is a normalization constant. By 
direct substitution we see that A annihilates ^o^) ^ annihilates ipo, consequently they 
are solutions with eigenvalue zero 

H.^ = A^(AiPq) = (17a) 
H + ^ = AiA^t) = (176) 
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but, these wave functions are not normalizable, i.e square integrable, at the same 
time, actually they need not be normalizable at all, but in the case when one of them is 
normalizable the supersymmetry is said to be unbroken. If non of them is normalizable, 
then non of the partner Hamiltonians has a zero eigenvalue and the supersymmetry is 
broken [HEIE]. 

From now on we will consider the case where the supersymmetry is unbroken, i.e 
one of and ip^ is normalizable, say ipQ. In this case H_ has an eigenvalue zero. We 
saw earlier that if H- has an eigenstate with a positive eigenvalue then H + also has the 
same eigenvalue and vice versa. Since the lowest eigenvalue for H_ is zero, then the 
lowest eigenvalue of H + is the first positive eigenvalue of H_. Moreover, the operators 
A and connect the eigenstates of H_ and H + except for the ground state of H 
[HI El H\ ■ This can be summarized in the following equations 

H_ip- = E-ip- n = 0,1,2,... (18a) 

H+i/>+ = E+i{>+ n = 0,1,2,... (186) 

E n+i = E n n = 0, 1, 2, and E~ = (18c) 

= —jL=A<l>Z+i n = 0,l,2,... (18d) 

E n+1 

^- +1 = -^=AV+ n= 0,1,2,... (18c) 



We notice here that the operators A and A^ annihilate and create nodes in wave 
functions, respectively, and this is a result of unbroken supersymmetry. They are not 
ladder operators though, because ladder operators move up and down the states of the 
same Hilbert space in contrast to these operators which shift us between different Hilbert 
spaces, for this reason some authors call them shift operators [TO"] . 

All the development so far depended on having a superpotential w(x) which gen- 
erates all the previous results. But the unbroken supersymmetric structure which has 
been considered so far can be constructed for any one dimensional quantum system with 
at least one bound state, i.e 

Every one-dimensional potential with bound states admits SUSY: 

Given a one dimensional quantum system with at least one bound state one can find 
a partner Hamiltonian which has exactly the same discrete spectrum except for the 
ground state energy of H_ [TTJ [i~2] . 

Let us consider having a one dimensional potential V- which admits bound state 
wave functions ip~, and the ground state energy is Eq , its Hamiltonian H_ given by 

H- = --^ + V4z)-EZ (19) 
can be factorized and be written as 

H = A ] A (20) 
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where A and are defined as in f 1 1 3 aj) and (113&P , respectively. The superpotential 

is 

w{x) = -— lnz/^ (21) 

where is the ground state wave function. The ground state energy of H_ is zero, 
and we can construct the supersymmetric partner Hamiltonian H + 

H + = AA^ (22) 

All the results of (18) are also true for our constructed partner Hamiltonians. 
This show us how to construct a supersymmetric system from a one dimensional 
potential. In fact this technique can be repeated again and again to obtain a sequence of 
Hamiltonians where all the adjacent pairs of Hamiltonians are supersymmetric partners, 
the cardinality of this sequence is the same as the dimensionality of the Hilbert space 
of H_, which is the number of bound states it possess [51 [TTJ, [13]. 



2. Shape Invariance 

We saw that supersymmetry allows us to construct and solve a hierarchy of Hamiltonians 
provided that we know the solution to one of the Hamiltonians in the sequence. 
The question now is, Can supersymmetry help us solve a supersymmetric system? 
The answer to this question was provided in 1983 by Gendenshtein [T3]. Having 
supersymmetry is not a sufficient condition for the system to be exactly solvable, 
because as we saw we can construct a supersymmetric Hamiltonian for any potential with 
bound state(s). Gendenshtein discovered another symmetry which if the supersymmetric 
system satisfies it will be an exactly solvable system, this symmetry is known as shape 
invariance. If our potential satisfies shape invariance we can readily write down its 
bound state spectrum, and with the help of the charge operators we can find the bound 
state wave functions pD, [TU, [15] . It turned out that all the potentials which were known 
to be exactly solvable until then have the shape invariance symmetry. 

If the supersymmetric partner potentials have the same dependence on x but differ 
in a parameter, in such a way that they are related to each other by a change of of that 
parameter, then they are said to be shape invariant. Gendenshtein stated this condition 
in this way, 

) = V-(x,a 1 ) + R(a 1 ) (23) 

where do is a parameter in our original potential whose ground state energy is zero. 
a\ = /(ao) where / is assumed to be an arbitrary function for the time being, and 
a n — f n {a-o) where f n is the composition of / with itself n times. The remainder R(ai) 
can be dependent on the parametrization variable do but never on x. In this case V_ is 
said to be shape invariant, and we can readily find its spectrum, take a look at H n and 
H n+ \ 

H n = -— + V4x,a n ) + J2 R ( a k) (24a) 
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, 2 n+1 

H n+1 = -— + V-(x,a n+1 ) + J2 R ( a k) (246) 

k=i 

d 2 n 
= --^ + V+(x,a n ) +J2 R ( a k) (24c) 

k=i 

where in the last step we applied ( |23i) . Here we see that H n and H n+ are 
super symmetric, so they have identical spectrum except for the ground state energy 
of H n which is, 

n 

EZ = J2R(*k) 

k=l 

but from the supersymmetry arguments discussed in section ([I]) we know that H n 
has the same spectrum of if _ = H° except for the first n levels of if_ which are missing 
from H n , so the ground state energy of the nth Hamiltonian H n is equal to the nth 
energy E~ of our original Hamiltonian So, the spectrum of H is 

n 

E~ = J2 R M E °=° ( 25 ) 

fc=l 

Supersymmetry together with shape invariance give us a formal expression for the 
wave functions ip~ of H_, for if we know the ground state wave function i/}q(x, ao) 
of H-(x,o,q), then the ground state wave function of H n (x,a n ) is ^(x,a n ), but the 
operator A*(x, a n -i) move us from the Hilbert space of H n to the Hilbert space of 
if n_1 , thus we can apply this operator repeatitively till we reach the Hilbert space of 
H = H°, we need to remember here that creates a node in the wave function each 
time we apply it, so applying it n times we will reach ip~(x, ao), hence 

ip~(x,a ) = A j! (x,a )A^(x,a 1 ) . . . A ] (x,a n _ 2 )A ] (x,a n ^ l )^(x 1 a n ) (26) 

where A^(x,a n ) is defined as 

A ] (x,a n ) = — + w(x,a n ) = 4~ - -j-lnipo(x,a n ) (27) 

{JjJL/ \JbtXj \JbJu 

This formal expression for the wave functions of shape invariant potentials were 
first presented in [16], the explicit expressions for the wave functions of all the well 
known shape invariant wave potentials have been worked out in [TTJ . Supersymmetry 
together with shape invariance can also be exploited to obtain the scattering matrices 



Different types of shape invariance have been studied in the literature, they are 
classified depending on the relation among the parameters which creates the shape 
invariant sequence of potentials. The first and the most dominant type is known as 
translational shape invariance, where the potential parameter is shifted, a\ = /(ao) = 
ao + a. This is the most well studied type of shape invariance, and many of the 
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generalizations of the shape invariance concept have been worked specifically for this 
type. In fact it turned out that all of the textbook examples of exactly solvable 
nonrelativistic one dimensional shape invariant quantum systems have this type of shape 
invariance. Complete lists of translational shape invariant potentials together with their 
wave functions and scattering matrices have been prepared (7[ [151 El IE] ■ 

Translational shape invariance was the only known shape invariance type for a 
long time, and some authors hypothesized that it is a necessary condition for shape 
invariance [20]. In 1993, as new types of shape invariance emerged, it was clear that 
translational shape invariance is not the only type. The first of the new types to emerge 
is the scaling shape invariance [2H [22], a\ = qa Q , < q < 1. There are also two 
nonlinear transformations which have been introduced together with the previous one 
in [22], a\ = qa$, p E Z, < q < 1, and ai = qaQ /{I + pao), < q,p < 1, pa Q <C 1. The 
eigenvalues, eigenfunction functions and transmission coefficients have been obtained 
algebraically for potentials with such shape invariance [2TJ [22] . 

Cyclic shape invariance has been introduced in the following years [23], [23] - In this 
type the parameters repeat themselves after cycle of p elements, a p = a . 

In all the types, except for the translational shape invariance, the potentials are 
not obtained in closed form, i.e in terms of elementary functions, they are obtained in 
series form. 

Another important generalization of the ordinary shape invariance relation ( 1231) is 
to have shape invariance in multi-steps. 

Shape invariance has also been generalized in another direction. In 1987, in an 
attempt to classify shape invariant potentials, Cooper et al.|25j introduced the idea of 
translational shape invariance in an "n" arbitrary but finite number of parameters, they 
were not successful at finding a solution for such a system, they were even pessimistic 
about the existence of a solution. It was an open problem for more than a decade, later 
on, by skillfully exploiting some properties of the Ricatti equation a solution was shown 
to exist and worked out[26j. 

Many successful efforts have been put to show the underlying algebraic structure 
of the shape invariance symmetry, the associated Lie algebras have also been identified 
[27] , and in the case of the non translational shape invariance nonlinear generalizations 
of Lie algebras have been obtained [2HI [29] . 

The technique of factorizing the Hamiltonian and using shape invariance to solve 
the schrodinger equation is not at all new, Schrodinger himself used this technique to 
solve the hydrogen atom, Dirac also used it to solve the harmonic oscillator. Of course 
the ideas of supersymmetry and the supergroups were not known back then. 

Later on, the technique of Schrodinger and Dirac was rejuvenated by Infeld and 
Hull and later they summarized their work in their infamous review |30j were they gave 
the name for the method, "The Factorization Method". They summarized the proce- 
dure of this technique to solve second order differential equations, they studied its range 
of applicability and classified the factorization types into six types (transformation be- 
tween the types exist), they also prepared a table for these six possible factorizations 
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which can be used to solve the differential equations simply by identifying to which 
factorization type it belongs [30]. For the historical progress of this method see [31] . 

Every Infeld and Hull factorizable potential is shape invariant but the con- 
verse is not true: 

As far as solving differential equations is concerned, the factorization technique offers 
almost the same power as the supersymmetry with translational shape invariance ap- 
proach. The supersymmetry and shape invariance offer physical insight and more group 
theoretic ideas of why such systems are exactly solvable. In fact shape invariance, in 
general, offers more than the factorization method, since the factorization method treats 
only the translational shape invariance. The equivalence of the factorization method to 
supersymmetry with translational shape invariance is easy to see, equation (3.1.2) in 
[30] is exactly equivalent to ( 1231) with translation of the shape invariance parameter 
[El [32]. 

The table prepared by Infeld and Hull, though complete for most purposes, is not 
the most general table possible, as they did not consider the most general solution of 
Ricatti equations. The most general solution is worked out in [32] . 

Shape invariance is a sufficient but not a necessary condition for exact solu- 
bility: 

Gendenshtein suggested that all exactly solvable potentials must be shape invariant 
[2] , but now the relation among exact solubility, shape invariance and supersymmetry 
is much clearer. One counter example is enough to refute Gendenshtein suggestions, 
and many such examples have been constructed [HI [25], [33], the most important of 
these counter examples is Natanzon class of potentials which are, in general, not shape 
invariant [25], [33] . 

3. Conclusion 

In this letter we defined, supersymmetric quantum mechanics, and shape invarinace. We 
also showed the relation which connects them to exact solubility and the factorization 
method. FigJT] summarizes these relations. 
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